Varying Alpha: New Constraints from Seasonal Variations 
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We analyse the constraints obtained from new atomic clock data on the possible time variation of 
the fine structure 'constant' and the electron- proton mass ratio and show how they are strengthened 
when the seasonal variation of Sun's gravitational field at the Earth's surface is taken into account. 
We compare these bounds with those obtainable from tests of the Weak Equivalence Principle and 
high-redshift observations of quasar absorption spectra consistent with time variations in the fine 
structure constant. 
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General relativity and the standard model of parti- 
cle physics depend on at least 27 independent parame- 
ters which determine the relative strengths of the differ- 
ent forces, matrix angles and phases, and the masses of 
all known fundamental particles. These parameters are 
commonly referred to as the fundamental constants of 
Nature and we are not able to explain or predict any of 
their precise numerical values. This characterisation may 
ultimately turn out to be a misnomer, as several modern 
proposals for the extension of fundamental physics be- 
yond the standard model predict that many of these pa- 
rameters are neither strictly fundamental nor constant. 
The precise values that the 'constants' take at a given 
instance are often determined in terms of the vacuum 
expectation values of one or more scalar fields and space- 
time variation of the 'constants', at some level, is there- 
fore a prediction common to most of these theories of new 
physics. We also recognize that if the true fundamen- 
tal theory exists in more than four spacetime dimensions 
then the most fundamental constants are defined there 
and the four-dimensional 'shadows' that we observe and 
call 'constants of Nature' can vary without undermining 
the status of constants in the higher-dimensional theory. 
Moreover, slow changes in the mean size of these extra 
dimensions will cause dimensionless constants to change 
at the same rate [lj. Consequently, observational inves- 
tigations of the constancy of these traditional constants 
may play an important role in elucidating the properties 
of any 'Theory of Everything'. 

Experimental and observational efforts to constrain the 
level of any possible time variation in some of these con- 
stants have a history that pre-dates the modern theo- 
ries which predict how they might vary (for overviews 
see Refs. [2]). Until recently, all observational studies 
found no evidence for variations in non-gravitational con- 
stants. In the last decade, however, data from a number 
of astronomical observations have provided suggestions 
that at least two of these constants, the fine structure 
constant: a — e 2 / 4ireoHc, and the electron-proton mass 
ratio: /i = m e /m p might have varied slightly over cos- 
mological time. Using a data set of 128 Keck-HIRES 
quasar absorption systems at redshifts 0.5 < z < 3, 



and a new many-multiplet (MM) analysis of the line 
separations between many pairs of atomic species pos- 
sessing relativistic corrections to their fine structure, 
Webb et al. Q found the observed absorption spectra 
to be consistent with a shift in the fine structure con- 
stant, a, between those redshifts and the present day, of 
Aa/a = a(z) - a(0)/a(0) = -0.57 ± 0.10 x 10" 5 . A 
smaller study of 23 VLT-UVES absorption systems be- 
tween 0.4 < z < 2.3 by Chand et al. [4( initially found 
Aa/a = — 0.6 ± 0.6 x 10~ 6 by using an approximate ver- 
sion of the full MM technique. However a recent reanal- 
ysis of the same data by Murphy et al. using the full un- 
biased MM method increased the uncertainties and sug- 
gested the revised figure of Aa/a — —0.64 ± 0.36 x 10~ 5 
for the same data [5j]. These investigations rely on the 
statistical gain from large samples of quasar absorption 
spectra. By contrast, probes of the electron-proton mass 
ration can use single objects effectively. Reinhold et al. 
[6J found a 3.5cr indication of a variation in the electron- 
proton mass ratio /i = m e /m pr o\ev the last 12Gyrs: 
A\i/\i = (—24.4 ± 5.9) x 10~ 6 from H 2 absorption in 
a different object at z = 2.8. However, recently Murphy 
et al 0] have exploited the /i sensitivity of ammonia in- 
version transitions [|| compared to rotational transitions 
of CO, HCN, and HCO + in the direction of the quasar 
B0218-I-357 at z = 0.68466 to yield a result that is con- 
sistent with no variation in fj, when systematic errors are 
more fully accounted for: Afi/fi = (+0.74 ± 0.47 sta t ± 
0.76 sys tcm) x 10~ 6 , corresponding to a time variation of 
A/m= (-l^iO.Sstatil^ystcmOxlO-^yr- 1 in the best 
fit ACDM cosmology. 

Any variation of a and /i today can also be constrained 
by direct laboratory searches. These are performed by 
comparing clocks based on different atomic frequency 
standards over a period of months or years. Until very 
recently, the most stringent constraints on the tempo- 
ral variation in a arose by combining measurements of 
the frequencies of Sr [|, Hg+ [nj, Yb+ [IH, and H 
relative to Caesium: a/a = (-3.3 ± 3.0) x lO^yr" 1 . 
Cingoz et al. also recently reported a less stringent 
limit of a/a = -(2.7 ± 2.6) x 10~ 15 yr _1 [15]; how- 
ever, if the systematics can be fully understood, an ul- 
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timate sensitivity of lO^^yr -1 is possible with their 
method [H)]. If a linear variation in a is assumed then 
the Murphy et. al. quasar measurements equate to 
a/a = (6.4 ± 1.4) x lO^yr" 1 3]. If the variation is 
due to a light scalar field described by a theory like that 
of Bekenstein [l3| and Sandvik, Barrow and Magueijo 
(BSBM) then the rate of change in the constants 
is exponentially damped during the recent dark-energy- 
dominated era of accelerated expansion, and one typi- 
cally predicts a/a = 1.1 ± 0.3 x 10~ 16 yr _1 from the 
Murphy et al data, which is not ruled out by the atomic 
clock constraints mentioned above. For comparison, the 
Oklo natural reactor constraints, which reflect the need 



for the Sm 1J ' J +?i — > Sm lj0 + 7 neutron capture resonance 
at 97.3 meV to have been present 1.8 — 2Gyr (z = 0.15) 
ago, as first pointed out by Shlyakhter [lj|, are currently 
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[17j Aa/a = (-0.8 ± 1.0) x 10~ s or (8.8 ± 0.7) x 10~ 8 
(because of the double-valued character of the neutron 
capture cross-section with reactor temperature) and [lit ] 
Aa/a > 4.5 x 10 -8 (6<r) when the non-thermal neutron 
spectrum is taken into account. However, there remain 
significant environmental uncertainties regarding the re- 
actor's early history and the deductions of bounds on 
constants. The quoted Oklo constraints on a apply only 
when all other constants are held to be fixed. If the quark 
masses to vary relative to the QCD scale, the ability of 
Oklo to constrain variations in a is greatly reduced [19j . 

Recently, Rosenband et al. [2(| measured the ratio 
of aluminium and mercury single-ion optical clock fre- 
quencies, /Ai+//Hg+, repeated over a period of about a 
year. From these measurements, the linear rate of change 
in this ratio was found to be (—5.3 ± 7.9) x 10 _17 yr _1 . 
These measurements provides the strongest limit yet on 
any temporal drift in the value of a: a/a = (—1.6 ± 



2.3) x 10 yr . This limit is strong enough to strongly 
rule out theoretical explanations of the change in a re- 
ported by Webb et al. in terms of the slow variation 
of an effectively massless scalar field, even allowing for 
the damping by cosmological acceleration, unless there 
is a significant new physical effect that slows the locally 
observed effects of changing a on cosmological scales (for 
a detailed analysis of global-local coupling of variations 
in constants see Refs. |2l|). 

It has been noted that if the 'constants' such as a or 
\x can vary, then in addition to a slow temporal drift 
one would also expect to see an annual modulation in 
their values. In many theories, the Sun perturbs the 
values of the constants by a factor roughly proportional 
to the Sun's Newtonian gravitational potential [22| (the 
contribution from the Earth's gravitational potential is 
about 14 times smaller than that of the Sun's at the 
Earth's surface). Hence the 'constants' depend on the 
distance from the Sun. Since the Earth's orbit around the 
Sun has a small ellipticity, the distance, r, between the 
Earth and Sun fluctuates annually, reaching a maximum 
at aphelion around the beginning of July and a minimum 



at perihelion in early January. It was shown in Ref. [23j] 
that in many varying constant models, the values of the 
constants measured here on Earth, would oscillate in a 
similar seasonal manner. Moreover, in many cases, this 
seasonal fluctuation is predicted to dominate over any 
linear temporal drift [231 ] . 

Specifically, let us suppose that the Sun creates a 
distance-dependent perturbation to the measured value 
of a coupling constant, C, of amplitude S In C = C{r). If 
this coupling constant is measured on the surface of an- 
other body (e.g. the Earth) which orbits the first body 
along an elliptical path with semi-major axis a, period 
T p , and eccentricity e <C 1, then to leading order in e, 



the annual fluctuation in C, SC, 



SC. 



annual 



c 



annual 
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will be given by 



= -c c cos( — ) +0(e 2 ), 



(1) 



where cq = eaC'(a), C (a) = dC (r) / dr\ r=a and t — nT p , 
for any integer n, corresponds to the moment of closest 
approach (perihelion). In the case of the Earth mov- 
ing around the Sun, over a period of 6 months from 
perihelion to aphelion one would therefore measure a 
change in the constant C equal to 2cc- Rosenband et 
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FIG. 1: (colour online). Frequency ratio fA\+/fn s + as mea- 
sured by Rosenband et al. [20]. The solid black line shows 
the maximum likelihood fit for a seasonal variation. 



al. [201 ] fitted a linear drift in a to their data finding 
a/a = (-1.6 ± 2.3) x 10~ 17 yr -1 . We fit the expected 
form of any annual fluctuation, Eq. [TJ to the measured 
values of /Ai+//Hg+- The data, taken from Ref. [20| . 
and our best fit line, are shown in Fig. [TJ It should 
be noted that the magnitude of the systematic errors for 
the middle three data points were not verified to the same 
jrecision as they were from the first and last data points 
20| . We do not expect to this have a great effect on 
the magnitude of the resulting constraint on k a . Using 
51n(/ A i+//H g +) = (3.19 + 0.008)fo/a, [U, a maximum 
likelihood fit to the data gives 



eaSa'(a) = (-0.89 ±0.84) x 10 
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(2) 
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where a = 149,597,887.5 km is the semi-major axis 
of the Earth's orbit, and Sa(r) is the perturbation in 
a due to the Sun's gravitational field. Assume that 
over solar system scales, the values of the scalar fields 
on which values of the 'constants' depend, vary with 
the local gravitational potential [22| . Hence, we have 
Sa(r)/a — k a AUQ(r), where k a is a theory-dependent 
multiplier, AUq is the change in the gravitational poten- 
tial of the Sun: J7©(r) = -GM /r, and so eaAU^a) — 
eGM Q /a = 1.65 x 1CT 10 . Hence, we find: 



k a = (-5.4 ±5.1) x 10" 



(3) 



The frequency shifts measured by Rosenband et al. [201 ] 
were not sensitive to changes in the electron-proton mass 
ratio: fi = m c /m p . Measurements of optical transition 
frequencies relative to Cs, Refs. 0, [l(l U, l3. are sensi- 
tive to both ji and a. H-maser atomic clocks [25j are also 
sensitive to variations in the light quark to proton mass 



ratio: q 



j/m p . We can use all these observations 



if we define two more gravitational coupling multipliers, 



ku, and k a 



and 5q/q = fc q A[/©. 



Refs. (a [lfl Hi] give k a + 0.36fc M = (-2.1 ± 3.2) x 10~ 6 , 
k a + 0.17fc M = (3.5 ± 6.0) x 10~ 7 , and k a + 0.13fc q = 
(1 ± 17) x 10~ 7 respectively. We also performed a boot- 
strap seasonal fluctuation fit (with 10 5 resamplings) to 



the Yb + frequency measurements of Peik et al. [ill [24 1 
giving fc Q + 0.51fc p = (7.1 ±3.4) x 10~ 6 . Combining these 
bounds with Eq. Q gives: 



fc M = (3.9 ±3.1) x 10~ fa , 
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(4) 

fc q = (0.1 ± 1.4) x 10~ b . (5) 

Recently, Blatt et al. combined data from measure- 
ments of H-maser [ii] and optical atomic clocks 0, [l(| , 
to bound the multipliers, k a , k^ and k q , finding: 

k a = (2.5 ±3.1) x 10~ 6 , 
fc M = (-1.3 ± 1.7) x 10" 5 , 
fc q = (-1.9 ±2.7) x 10~ 5 . 

The constraint on k„ derived in this paper from the data 
of Rosenband et al. [201 ] therefore represents an improve- 
ment by about two orders of magnitude over the previous 
best bound. This improved bound on fc Q combined with 
data found by Peik et al. [ll|, |2J] has also produced an 



order of magnitude improvement in the determination of 
fc M and a slight improvement in the constraint on k q . 

Seasonal fluctuations are predicted by a varying con- 
stant theory because the scalar field which drives the 
variation in the constant couples to normal matter. The 
presence of the Sun therefore induces gradients in scalar 
field and the 'constants', and it is essentially these gradi- 
ents that are detectable as seasonal variables. However, 
as is well known, gradients in a scalar field which cou- 
ples to normal matter result in new or 'fifth' forces with 
pseudo-gravitational effects. In the case of varying a and 



fj, theories, these forces are almost always composition de- 
pendent, which would violate the universality of free-fall 
and hence the weak equivalence principle (WEP). The 
magnitude of any composition-dependent fifth force to- 
ward the Sun is currently constrained to be no stronger 
than 10~ 12 — 10 -13 times than the gravitational force 
[26j]. In the context of a given theory the constraints 
from WEP tests indirectly bound k a . Indeed, they often 
provide the tightest constraints on fc Q 0, 0, [5§] . 

A recent and thorough analysis of the WEP violation 
constraints on k a 28|, SO] found: 



k a = (0.3 ±1.7) x 10" 9 , 

with a similar constraint on k q . It must be noted, how- 
ever, that this result is still subject to theoretical un- 
certainty, especially regarding the dependence of nuclear 
properties on quark masses. For instance, it was also 
noted in Ref. [28| that if certain (fairly reasonable) as- 
sumptions about nuclear structure are dropped, the la 
error bars on k a increase by about an order of magni- 
tude to: ±1.4 x 10~ 8 . Despite these uncertainties, for 
many theories of varying a, WEP violation constraints 
from laboratory experiments or lunar laser ranging [27| 
still provide the strongest, albeit indirect, bound on fc Q . 

We have shown above that direct constraints on any 
change in a that is proportional to the local gravitational 
potential (as in general theories of its spacetime variation 
[22j) are now within less than 1.5 orders of magnitude of 
those extrapolated, under certain fairly reasonable as- 
sumptions, from the non-detection of any WEP viola- 
tion. One might wonder how much further the sensitiv- 
ity of atomic clock experiments to variations in a would 
have to improve before direct constraints on k a would 
surpass those coming from the current WEP violation 
bounds. Suppose that over a few days, the ratio of two 
transition frequencies, /a//b, can be measured with a 
sensitivity at, and that <51n(/A//B) = S a Sa/a (typically 
S a ~ 0(1), although some transitions exhibit a greatly 
increased sensitivity [3ll]). The sensitivity to changes in 
a is then given by a a — ai/ S a . By simulating data sets, 
we found that the sensitivity to fc Q is significantly im- 
proved if one makes N m ^12 measurements per year (at 
roughly regular intervals). With 7V m > 12, by performing 
a bootstrap linear regression with 10 5 resamplings of the 
simulated data points, we find that the sensitivity, a^ to 
k a is roughly: 

cr k w 0.69 x 10 r 



y/N y (N m - 1) 

where N y is the number of years for which data is 
taken. The total number of measurements is therefore 
N y N m . Indirect constraints currently have a sensitivity 
ofcik = 1.7xl0 -9 [il]. This would be surpassed by direct 
measurements if a a < 2.5y/N y (N m — 1) x 10~ 19 . For ex- 
ample, if we take measurements every 20 days (or so) over 
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a single year (N m — 18, N y = 1) then we would require 
a a < 1CP 18 . The experiment conducted by Rosenband 
et al. [20| currently has a <T[ « 5 x 10~ 17 and S a = 3.2, 
so a a ~ 1.6 x 10 -17 ; one would therefore require an un- 
realistic number of measurements N y (N m — 1) > 4100 
to achieve the desired sensitivity to k a . Cingoz et al. 
obtained constraints on k a by monitoring the transition 
frequencies between nearly degenerate, opposite-parity 
levels in two isotopes of atomic Dysprosium [HI]. Cur- 
rently sensitively to frequency changes is 7.5Hz, giving 
(T a = 5 x 10~ 15 . However, an ultimate sensitivity of 
about a mHz is predicted to be achievable [32^ . This 
would give a a w 7 x 10~ 19 and one would only require 
N y (N m - 1) > 8 to surpass the current WEP viola- 
tion constraints on k a . In addition, Flambaum noted 
in Ref. [3~0 ] that the transition between the ground and 
the first excited states in the 229 Th nucleus is particu- 
larly sensitive to changes in a and fi with S a ~ 10 5 . 
By measuring this transition frequency relative to, say, 
a Cs frequency standard, and without any improvement 
in measurement accuracy over what has currently been 
achieved, it may be possible bring a a down to about 
10~ 20 , giving cr k < 2 x 10~ n with N m > 12 - two orders 
of magnitude better than the current WEP violation con- 
straints. The motivation for a future space-based test of 
the WEP with possible sensitivity of order 10~ 15 — 10~ 18 
therefore remains very great 33j, [34 1 . 

In summary: we have shown how new laboratory con- 
straints on possible time variation in the fine structure 
'constant' and the electron-proton mass ratio can yield 
more sensitive limits by incorporating the effects of the 
seasonal variation of the Sun's gravitational field at the 
Earth's surface. This seasonal variation is expected in all 
theories which require that the covariant d'Alembertian 
of any scalar field driving variation of a 'constant' is pro- 
portional to the dominant local source of gravitational 
potential 29]. The recent experimental results from 
Rosenband et al [2(| and Peik et al, 11 1 have reached 
the sensitivity of the quasar observations of varying a 
and /j, made at high redshift and we have shown may 
soon provide stronger bounds on varying constants than 
conventional ground-based WEP experiments. 
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